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ABSTRACT 

The clustering of X-ray selected AGN appears to be a valuable tool for extracting cosmological in- 
formation. Using the recent high-precision angular clustering results of ^ 30000 XMM-Newton soft 
(0.5-2 keV) X-ray sources (Ebrero et al.), which have a median redshift of z ~ 1, and assuming a flat 
geometry, a constant in comoving coordinates AGN clustering evolution and the AGN bias evolution 
model of Basilakos et al., we manage to break the ftm — crs degeneracy. The resulting cosmological 
constraints are: 0^, = O-'^'^-o 'ot^ "O-^Olg^ig and as = ^■'^^-o.it' while the dark matter host halo 
mass, in which the X-ray selected AGN are presumed to reside, is M = 2.50l°'5g ^ W^^h~^MQ. For the 
constant A model (w=— 1) we find fi^ = 0.24 ± 0.06 and as — 0.83lg ijg, in good agreement with recent 
studies based on cluster abundances, weak lensing and the CMB, but in disagreement with the recent 
bulk flow analysis. 
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1. INTRODUCTION 

A large variety of cosmologically relevant data, based 
on the combination of galaxy clustering, the supernova la's 
Hubble relation, the cosmic microwave background (CMB) 
fluctuations and weak-lensing strongly support a flat uni- 
verse, containing cold dark matter (CDM) and "dark en- 
ergy" which is necessary to explain the observed acceler- 
ated cosmic expansion (eg., Komatsu et al. 2010; Hicken 
et al. 2009; Fu et al. 2008 and references therein). 

The nature of the mechanism that is responsible for the 
late-time acceleration of the Hubble expansion is a fun- 
damental problem in modern theoretical physics and cos- 
mology. Due to the absence of a physically well-motivated 
fundamental theory, various proposals have been suggested 
in the literature, among which a cosmological constant, 
a time varying vacuum quintessence, fc— essence, vector 
fields, phantom, tachyons, Chaplygin gas, etc (eg., Wein- 
berg 1989; Peebles & Ratra 2003; Boehmer & Harko 2007; 
Padmanabhan 2008 and references therein). Note that 
the simplest pattern of dark energy corresponds to a scalar 
field having a self-interaction potential with the associated 
field energy density decreasing with a slower rate than the 
matter energy density. In such case the dark energy com- 
ponent is described by an equation of state pg — wpQ with 
w< —1/3 (dubbed "quintessence", eg. Peebles & Ratra 
2003 and references therein). The traditional cosmologi- 
cal constant (A) model corresponds to w= —1. The via- 
bility of the different dark-energy models in reproducing 
the current excellent cosmological data and the require- 
ments of galaxy formation is a subject of intense work (eg. 
Basilakos, Plionis & Sola 2009 and references therein). 

Another important cosmological parameter is the nor- 
malization of the cold dark matter power spectrum in the 
form of the rms density fluctuations in spheres of radius 
8h^^ Mpc, the so-called erg. There is a degenerate relation 



between as and (eg. Eke, Cole & Frenk 1996; Wang & 
Steinhardt 1998; Henry et al. 2009; Rozo et al. 2009 and 
references therein) and it is important to improve current 
constraints in order to break such degeneracies. Further- 
more, there are also apparent inconsistencies between the 
values of as provided by different observational methods, 
among which the most deviant and problematic for the 
concordance cosmology, is provided by the recent bulk flow 
analysis of Watkins, Feldman & Hudson (2009). 

In this paper we extend our previous work (Basilakos & 
Plionis 2009; hereafter BP09), using the angular clustering 
of the largest sample of high-z X-ray selected active galac- 
tic nuclei (Ebrero et al. 2009a), in an attempt to break 
the as — degeneracy within spatially flat cosmological 
models. 

2. BASIC METHODOLOGY 

The main ingredients of the method used to put cosmo- 
logical constraints based on the angular clustering of some 
extragalactic mass-tracer, has been already presented in 
our previous papers (see also Matsubara 2004; BP09 and 
references therein). It consists in comparing the observed 
angular clustering with that predicted by different primor- 
dial fluctuations power-spectra, using Limber's equation 
to invert from spatial to angular clustering. By minimiz- 
ing the differences of the observed and predicted angular 
correlation function, one can constrain the cosmological 
parameters that enter in the power-spectrum determina- 
tion as well as in Limber's inversion. Below we present 
only the main steps of the procedure. 

2.1. Theoretical Angular and Spatial Clustering 

Using the well known Limber's inversion equation (Lim- 
ber 1953), we can relate the angular and spatial clustering 
of any extragalactic population under the assumption of 
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power-law correlations and the small angle approximation 
(see details in BP09). After some algebraic calculations 
and within the context of flat spatial geometry, we can 
easily write the angular correlation function as: 

w{e)^2^j^ (^-— j E{z)dzj^ e(r,z)du , (1) 

where AN /dz is the source redshift distribution, estimated 
by integrating the appropriate source luminosity function 
(in our case that of Ebrero et al. 2009b), folding in also 
the area curve of the survey. We also have 

E{Z) = [nmil + Zf + (1 - + z)3(l+-)]l/2 ^ (2) 

with w the dark-energy equation of state parameter given 
by pQ = wpQ with w< —1/3. The source spatial correla- 
tion function is: 

f(r,z) = (l + z)-(3+^)62(z)^DM(r), (3) 
where b{z) is the evolution of the linear bias factor, e 
is a parameter related to the model of AGN clustering 
evolution (eg. de Zotti et al. 1990 j3 and Cdm(?') is the 
corresponding correlation function of the underlying dark 
matter distribution, given by the Fourier transform of the 
spatial power spectrum P{k) of the matter fluctuations, 
linearly extrapolated to the present epoch: 



-dk 



(4) 



JO kr 
We use the nominal functional form of the CDM power 
spectrum, P{k) = Pok'^T'^ik), with T{k) the CDM trans- 
fer function (Bardeen et al. 1986; Sugiyama 1995) and 
n ~ 0.96, following the 5 (and 7)-year WMAP results 
(Komatsu et al. 2010), and a baryonic density of fli^h^ = 
0.022(±0.002). The normalization of the power-spectrum, 
Pq, can be parametrized by the rms mass fluctuations on 
i?8 = Sh^^Mpc scales (cts), according to: 

-1 



^^0 



Zn (7g 



T^{k)k''+'W^{kRs)dk 



(5) 



where W{kRs) = 3(sinfci?8 - kRsCOskRs)/ikRs)^. Re- 
garding the Hubble constant we use either Hq — 71 
kms-^Mpc-i (Freedman 2001; Komatsu et al. 2010) or 
Ho ~ 74 kms-^Mpc-i (Riess et al. 2009). Note, that in 
the current analysis we also utilize the non-linear correc- 
tions introduced by Peacock & Dodds (1994). 

2.2. X-ray AGN bias evolution 

The notion of the bias between mass-tracers and un- 
derlying DM mass is an essential ingredient for CDM 
models in order to reproduce the observed extragalactic 
source distribution (eg. Kaiser 1984; Davis et al. 1985; 
Bardeen et al. 1986)Although a large number of models 
have been proposed in the literature to model the evolu- 
tion of the bias factor, in the current analysis we use our 
own approach, which was described initially in Basilakos 
& Plionis (2001; 2003) and extended in Basilakos, Plionis 
& Ragone-Figueroa (2008; hereafter BPR08). 

For a benefit of the reader we remind that our bias model 
is based on linear perturbation theory and the Friedmann- 
Lemaitre solutions of the cosmological field equations, 

* Following Kiindic (1997) and Basilakos & Plionis (2005; 2006) we 
^ For the benefit of the reader we present the values of the constants 
while in eq.ltSj we have A = 5 X 10~^ and m = 2.62 for 1 X 10^'^ < 
1 X 10" < M/h-^MQ < 6 X 10". 



while it also allows for interactions and merging of the 
mass tracers. Under the usual assumption that each X-ray 
AGN is hosted by a dark matter halo of the same mass, 
we can present analytically its bias evolution behavior. A 
more realistic view, however, of the AGN host halo hav- 
ing a spread of masses around a given value, with a given 
distribution that does not change significantly with red- 
shift, should not alter the predictions of our bias evolution 
model. 

For the case of a spatially flat cosmological model, our 
bias evolution model predicts: 

6(Af, z) = Ci{M)E{z) + C2{M)E{z)I{z) + y{z) -f 1 , (6) 

where 



y{z) = E{z) 



lC{x)I{x)dx 



I{z) 



lC{x)dx 



(7) 



(8) 
(9) 



{l + xf ' ' J„ (l + x)^ 

with JC{z) = f{z)E^{z), I{z) = /;°(1 + xfdx/E'^ix) 

f(z) = A{m-2){l + 7TE{z)/D{z) , 

Cia{M) ~ ai,2(M/lO"/i-iM0)''i.= , 

The various constants are given in BPRO^. Note that 
D{z) is the linear growth factor (scaled to unity at the 
present time), useful expressions of which can be found 
for the dark energy models in Silveira & Waga (1994) and 
in Basilakos (2003). 
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Fig. 1. — The observed evolution of AGN bias (different points) 
compared with the BPR08 model predictions (curves). Optically 
selected SDSS and 2dF quasars are represented by empty dots 
and crosses, respectively, while X-ray selected AGN by filled (blue) 
squares. In the insert we plot the most recent optical QSO bias 
values based on the SDSS quasar uniform sample (Ross et al. 
2009). The curves represent the expectations of the BPR08 model, 
with the solid (red) lines corresponding to a DM halo mass of 
M = 10" h-^ Mq and the dashed line to M = 2.5 X lO" Mq. 

In order to provide an insight on the success or failure 
of our bias evolution model, we compare in Fig.l the mea- 
sured bias values of optical and X-ray selected AGNs with 
our 6(z) model. The bias of optical quasars by Groom 
et al. (2005), Myers et al. (2007), Shen et al. (2007) 
and Ross et al. (2009) based on the 2dF QSOs (open 
circles), SDSS DR4 (crosses), SDSS DR5 (sohd points) 
and the SDSS quasar uniform sample (insert panel), are 

use the constant in comoving coordinates clustering model, ie., e = —1.2. 

for a few selected cases: ai = 3.29, /3i = 0.34 and 02 = —0.36, /32 = 0.32, 
Af/h-^M© < 3 X 10", and A = 6 X 10~^ and m = 2.54 for the case of 
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well approximated by our b{z) model for a DM halo of 
10^^h~^ Mq (solid red line) in agreement with previous 
studies (Porciani, Magliocchetti, Norberg 2004; Croom et 
al. 2005; Negrello, Magliocchetti & de Zotti 2006; Hop- 
kins et al. 2007). However, what is worth stressing is that 
our model is (to our knowledge) the only one that can si- 
multaneously fit the lower redshift {z < 2.5) optical AGN 
bias with the higher (z > 3) results of Shen et al. (2007) 
for the same halo mass of lO^^^h^-^ Mq. The solid (blue) 
squares represent the bias of the soft X-ray selected AGNs, 
based on a variety of X-ray surveys (eg. Basilakos et al. 
2005; Puccetti et al. 2006; Gilli et al. 2009; Ebrero et al. 
2009a). The model b{z) curve (dashed line) that fits these 
results correspond to halo masses M = 2.5 x 10^'^ Mq, 
strongly indicating that X-ray and optically selected AGN 
do not inhabit the same DM halos. 

3. COSMOLOGICAL PARAMETER ESTIMATION 

We use the most recent measurement of the angular cor- 
relation function of X-ray selected AGN (Ebrero et al. 
2009a). This measurement is based on a sample (here- 
after 2XMM) constructed from 1063 XMM-Newton obser- 
vations at high galactic latitudes and includes ~ 30000 
soft (0.5-2 keV) point sources within an effective area of 
~ 125.5 deg^ and an effective fiux-limit of fx > 1.4 x 10"^^ 
erg cm~^ s^^ (for more details see Mateos et al. 2008). 
Notice that the redshift selection function of the X-ray 
sources, obtained by using the soft-band luminosity func- 
tion of Ebrero et al. (2009b) , that takes into account the 
realistic luminosity dependent density evolution of the X- 
rays sources, predicts a characteristic depth of z ~ 1. 

In BP09, using the 2XMM clustering, we already pro- 
vided stringent cosmological constraints in the flm — w 
plane, using as a prior a fiat cosmology and the WMAP7 
power-spectrum normalization value of Komatsu et al. 
(2010). In the current analysis we relax the latter prior 
and allow cg to be a free parameter to be fitted by the 
data. Therefore the corresponding free-parameter vector 
that enters the standard likelihood procedure, which 
compares the observed and predicted clustering, is: p = 
(il,„, w, CTg, M), with M the AGN host dark matter halo 
mass, which enters in our BPR08 biasing evolution scheme. 

The likelihood estimatoi0, is defined as: >Cagn(p) oc 
exp[-xiGN(p)/2] with: 



Xagn 



(p) = J2 Kh(0.,p) - woMf ' (10) 



where n and Ci is the number of logarithmic bins (n = 13) 
and the uncertainty of the observed angular correlation 
function respectively, while ag- corresponds to the width 
of the angular separation bins. 

We sample the various parameters in a grid as follows: 
the matter density flm & [0.01, 1] in steps of 0.01; the equa- 
tion of state parameter wG [—1.6, —0.34] in steps of 0.01; 
the rms matter fiuctuations erg £ [0.4, 1.4] in steps of 0.01 
and the parent dark matter halo AZ/lO^^/i^^M© e [0.1,4] 
in steps of 0.1. Note that we have allowed the parameter 
w to take values below -1. 

Our main results are listed in Table 1, where we quote 
the best fit parameters with the corresponding la uncer- 
tainties, for two different values of the Hubble constant. 

^ Likelihoods are normalized to their maximum values. 



Small variations around ^ 71 kms^^Mpc"^ (which is the 
value used in the rest of the paper), appear to provide sta- 
tistically indistinguishable results. The likelihood func- 



7+0.03 



tion of the soft X-ray sources peaks at i7„j = 0.27_(j Q5, 



w= -0.90to,\l, cTs = 0.74+^:i^ and M = -^-^l^U x 
10^^ h~^MQ, with a reduced of ~ 4. Such a large x^/df 
value is caused by the measured small w{9) uncertainties 
in combination with the observed w{0) sinusoidal modu- 
lation (see BP09). Had we used a 2<t w{6) uncertainty 
in eq. ljlOp we would have obtained roughly the same con- 
straints and a reduced of ^ 1. The apparent sinusoidal 
w{9) modulation is a subject of further investigation. 
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Fig. 2. — Likelihood contours (Icr, 2cr and Scr) in the following 
planes: (Om,w) (upper left panel), {Qrn,crg) (upper right panel), 
(o-gjw) (bottom left panel) and {as,M) (bottom right panel). In 
the upper two panels we show for clarity our current solution with 
thick (red) contours, while the dashed contours correspond to our 
previous analysis, based on the shallower XMM/2dF survey (Basi- 
lakos & Plionis 2006). 

In Fig. 2 we present the Icr, 2a and 3cr confidence lev- 
els (corresponding to where — 21n£/£„iax equals 2.30, 6.16 
and 11.83) in the (17m,w), (f2m,crg), (CTg,w) and (crgjM), 
planes, by marginalizing the first one over M and erg, the 
second one over M and w, the third one over M and rim 
and the last over n,„ and w. We also present, with dashed 
lines, our previous solution of Basilakos & Plionis (2006), 
which where derived by using the shallower (effective fiux- 
limit of fx > 2.7 X 10"^'' erg cm^^ s^^) and significantly 
smaller (-- 2.3 deg^) XMM/2dF survey (Basilakos et al. 
2005). Comparing our current results with our previous 
analysis it becomes evident that with the current high- 
precision X-ray AGN correlation function of Ebrero et al. 
(2009a) we have achieved to break the ftm — erg degener- 
acy and to substantially improve the constraints on $7^, 
w and (Tg. However, there are still degeneracies, the most 
important of which is in the w— cg plane. 

It should be mentioned that some recent works, based on 
the large-scale bulk flows, strongly challenge the concor- 
dance ACDM cosmology by implying a very large CTg value. 
Indeed, Watkins et al. (2009), using a variety of tracers to 
measure the bulk flow on scales of ~ 100 h^^Mpc, found 
a value of ~ 400 km s~^ that implies a cg normalization 
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which is a factor of ^ 2 larger than what expected in 
the concordance cosmology. On the high erg side are also 
the results of Reichardt et al. (2009), based on the sec- 
ondary Sunayev-Zeldovich anisotropics in the CMB, pro- 
viding (Tg ~ 0.94, as well as a novel analysis based on the 
integrated Sachs- Wolfe effect (Ho et al. 2008). 

Contrary to the above results, our X-ray AGN cluster- 
ing analysis provides a ag value consistent with the con- 
cordance cosmology and in agreement with a variety of 
other studies. In particular, for w= —1 (A cosmology) 
and M = 2.5toi x lO^"'' h-^M^, we find 0„ = 0.24 ±0.06 
and ag = O.SS^q :[g (sec Table 1). Our results arc in 
agreement with those of recent cluster abundances stud- 
ies, providing (for w= -1): erg = 0.86 ± 0.04(1^^/0.32) 
(Henry et al. 2009) and <Tg = 0.83 ± 0.03(f2„/0.25)-0 '*i 
(Rozo et al. 2009). Furthermore, Mantz et al. (2009) 
using as a new cosmological tool the simultaneous fit of 
the cosmological parameter space and the cluster X-ray 
luminosity-mass relation, also broke the 12^ — ag degen- 
eracy and found: fi™ = 0.23 ± 0.04, trg = 0.82 ± 0.05 
and w= -1.01 ± 0.20 (sec their Table 2 and Fig. 4), 
which are in excellent agreement with our w= —1 results 
(see last row in our Table 1). Note, that their combined 
analysis (utilizing also the CMB, baryonic acoustic oscil- 
lations and gas mass fraction) provides = 0.27 ± 0.02, 
w= -0.96 ± 0.06 and erg = 0.79 ± 0.03. Moreover, Fu et 
al. (2008) based on a weak-lensing analysis found the de- 
generate combination tjg = 0.837 ± 0.084(O„/0.25)° -S3. 
From the peculiar velocities statistical analysis. Pike & 
Hudson (2005) and Cabre & Gaztahaga (2009) obtained 
o-g = 0.88 ± 0.05(fJ„/0.25)-°-53 and ag = 0.85 ± 0.06 
(for Q.m = 0.245), respectively. The consistency of all the 
(seven) previously mentioned works (including the current 



study) can be also appreciated from their average erg value 
which is (for w= —1 and ignoring the different val- 
ues): (erg) = 0.844 ± 0.009, where the quoted uncertainty 
is the Icr scatter of the mean. Note that the combined 
WMAP 7-years+SNIa+BAO analysis of Koniatsu et al. 
(2010) provide a slightly lower value of ag = 0.809 ±0.024 
(with 0„ = 0.272 ± 0.015 and w= -0.98 ± 0.05). 

4. CONCLUSIONS 

We have used the recent angular clustering measure- 
ments of high-z X-ray selected AGN, identified as soft 
(0.5-2 keV) XMM point sources (Ebrero et al. 2009a), 
in order to break the degeneracy between the rms mass 
fluctuations ag and Applying a standard likelihood 
procedure, assuming a constant in comoving coordinates 
AGN clustering evolution, the bias evolution model of 
Basilakos et al. (2008) and a spatially flat geometry, we 
put relatively stringent constraints on the main cosmologi- 
cal parameters, given by: D.„, = O.27to !5-^, w= -0.90to Jg 

and (Tg — 0.74^0 'ijj- We also find that the dark matter 
host halo mass, in which the X-ray selected AGN are as- 
sumed to reside, is M = 2.50l?:|g x 10^^/i"^Mq. Finally, 
if we marginalize over the previous host halo mass and 
w= —1 (A cosmology), we find flm = 0.24 ± 0.06 and 
erg = 0.83tQ;ig. 
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Table 1 

The best fit cosmological parameters from the likelihood analysis. 



iJo/kms-'Mpc"' 




w 


0-8 


M/lO^^/j-i^^ 


71 
74 


n 97+().(« 

"•^'-0.05 
U.ZD_o,05 


-0 9G+"-i" 

^••^"-0.16 

-0.92^^4 


^- '^-0.12 

0.72l°;l^ 


^•^^-1.50 

2 (^n+o-^o 

^•^^-1.50 


71 


0.24 ±0.06 


-1 


0.83t°;lJ 


2.50 



"Errors of the fitted parameters represent 1<t uncertainties. 



